Brownian motion meets Riemann curvature 



Pavel Castro VillarreaE 

Centra de Estudios en Ffsica y Matematicas Basicas y Aplicadas, 
Universidad Autonoma de Chiapas, C.P. 29000, 1428 Tuxtia Gutierrez, Chiapas, MEXICO 



O Abstract 

(N 

K^ The general covariance of the diffusion equation is exploited in order to explore the curvature 

C^ effects appearing on brownian motion over a d-dimensional curved manifold. We use the local frame 

^ defined by the so called Riemann normal coordinates to derive a general formula for the mean- 

square geodesic distance (MSD) at the short-time regime. This formula is written in terms of 0{d) 

invariants that depend on the Riemann curvature tensor. We study the n-dimensional sphere case 

'~~P to validate these results. We also show that the diffusion for positive constant curvature is slower 

Q than the diffusion in a plane space, while the diffusion for negative constant curvature turns out to 

^ be faster. Finally the two-dimensional case is emphasized, as it is relevant for the single particle 



a 



diffusion on biomembranes. 



-(— > 

^ 1 Introduction 

2 The brownian motion phenomenon occurs in a wide diversity of physical areas; from colloidal physics 

F^ to quantum gravity and biophysics (see for instance [lj,[2] and |3]). In the last decade, motivated by 

ti problems in biophysics (see [4l[5]), an intense activity has emerged on the study of diffusion processes on 

Q curved manifolds. For instance, the transport phenomena occurring in cell membranes is an interesting and 

'— ' complex problem. In particular, the random motion of a single integral protein or lipid in a cell membrane 

^H is difficult to realise, mainly because of the interactions and obstacles with the remainder components of 

^ the cell. In addition to this difficulty, the thermal fluctuations produce shape undulations on the curved 

ly-j membrane [51 [71 [B]. By simplifying this problem to the study of free diffusion on a cell membrane, 

^ considered as a regular and continuos two-dimensional surface, Smolouchovski's diffusion equation has 

• been proposed in [9, JOJ. Explicit formulas for constant mean and gaussian curvatures have been presented 

Q at [TT] and [12], respectively. Nevertheless, the analytical issues that appear on the general surface case, 

O have motivated the incorporation of novel computer simulations [121 [H] (see also [5] for a related work). 

T1 From a theoretical stand point, the brownian motion can be used to probe the geometry of the manifold 

.^ in the spirit of Kac's famous question: Can one hear the shape of a drum? [15j, and vice versa, the 

^ geometry will cause a change in the standard randomness of the particle motion [16J. It is then imminent 

^ to analyse a quantitative contribution coming from the geometry, and to be specific, how the curvature of 

the manifold affects the motion. These questions have been arisen in [12J, where the local concentration 

of a diffused substance is obtained in terms of the local curvature. Furthermore, a detailed study of the 

way on how the mean-values are affected by the curvature for two-dimensional manifolds is presented in 
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|17) . with a special emphasis on developable and isotropic surfaces. Also, the study of particular cases 
presented in [18J suggests that for negative Gaussian curvature, the diffusion accelerates; whereas, the 
diffusion decelerates for surfaces with positive Gaussian curvature. We should point out that, although 
with different formalism, these same questions have been posted in [I^ . 

In this work we explore the curvature effects on the brownian motion when the particle movement 
takes place on a d-dimensional riemannian manifold. These effects manifest differently for different physical 
observables. Here, we will use the geodesic distance as the displacement of the particle, but some other 
observables can also be defined [lOl [T31 [IT], where either intrinsic or extrinsic properties of the manifold 
can be probed [20] ■ Furthermore, we will take advantage of the general covariance of the diffusion 
equation to use a special frame defined by the so called Riemann Normal Coordinates (RNC) [2]. In this 
frame we will compute curvature corrections for the mean-square geodesic distance. In particular, we use 
a technique developed in [22J, originally used to compute curvature corrections that appear in effective 
actions of field theory on curved spaces (see for instance [23, 24j). Related work concerning the RNC can 
be found at [251 |2S]- It is remarkable that, using these coordinates, the geodesic curves on the manifold 
look like stright lines and therefore, the square geodesic distance s^ will have the same structure as the 
square distance in an euclidean space: s^ = 5aby""y'' (where y"- are the RNC). As it is shown in this paper, 
the mean-square geodesic distance is clearly isometric, this is because it only depends on 0(d) invariant 
combinations of the Riemann curvature tensor. 

This paper is organized as follows: In section 2 we summarize the geometrical concepts used to 
approach the brownian motion. In particular, we introduce the frame defined by the Riemann Normal 
Coordinates. In section 3, we present the diffusion equation on curved manifolds and we derive general 
remarks concerning the short time regime. In section 4, we focus on the computation of the mean-square 
geodesic displacement using RNC. In particular, we study the curvature effects for the manifolds with 
constant curvature and the two-dimensional case. The n-dimensional sphere is also explore to validate 
our results. Finally, in section 5 we sintetize the main conclusions and perspectives of this work. 

2 Geometrical preliminaries and notation 

In this section we review the preliminary notions about manifolds and riemannian geometries (following 
[2ZJ) needed to describe the brownian motion. Let us call M a d-dimensional manifold with local coordi- 
nates (/) : f/ C M — i- W^, where f/ is a local neighbourhood and MJ^ is the rf-dimensional euclidean space. 
As a consequence of the differentiablity of the map ip, the set U is locally diffeomorphic to a piece of the 
euclidean space. The local coordinates are also denoted by (^{p) = {x^, ■ ■ ■ ,x^) where p G M. For each 
point p on the manifold we associate a vector space called the tangent space TpM, whose elements are 
denoted by capital letters X,Y, Z,- ■ ■ . 

We are interested in manifolds endowed with a riemannian metric. \f gp : TpM x TpM — )■ M denotes 
a riemannian metric, we write 

gp = gab dx" (g) (ix^ (1) 

where {dx""} constitute a one-form basis of the dual tangent space T*M and gah is the metric tensor. 
The Einstein summation rule is adopted for the repeated indexes. The knowledge of the metric tensor 



components allow us to compute further geometrical quantities. The geometrical meaning of how the 
manifold is curved is defined in terms of the torsion tensor and the Riemann curvature tensor [57], 

T(X,r) = VxY-VyX-[X,Y], 
RiX,Y,Z) = VxVYZ-VYVxZ-Vyx,Y]Z (2) 

where V is the affine connection. For our purposes, we will use the Levi-Civita connection, which is the 
only one compatible with the metric. Using the coordinate basis {ca} = {da} of the tangent space TpM, 
the affine connection defines the components T\^ by 

VaCfe = Ve^Cfe = T^f^Cc, (3) 

where Va stands for the covariant derivative and for the Christoffel symbols r"j,^F| The torsion is a (1,2) 
type tensor T = T^^ e^ ® dx^ ® dx'' whereas the Riemann curvature is a (1,3) type one R = R\^ €„ ® 
dx^®dx'^®dx'^. The components of these tensors are given in terms of the Chrystoffel symbols^]. Clearly, 
the manifold is free torsion for the Levy-Civita connection, because the Christoffel symbols are symmetric. 
The quantity Rabcd = 9afR\cd satisfies the following useful identities Rabcd = -Rbacd = -Rabdc = Rcdab- 
Using the Riemann curvature components we can defined the Ricci tensor Rab = R'^acb ^""^ ^^e scalar 
curvature Rg = g'^^Rab- 

There is an important device introduced by Riemann, nowadays called the Riemann Normal Coordinates 
pi] . This coordinate system can be defined by mapping a point p on the manifold to the origin of W^ 
and the following conditions 

QabiS^) = Sab, y^9ab{y) = V'^^ab- (4) 

As it is point out at [25], the second condition is equivalent to the following gauge condition on the affine 
connection (see appendix A) 

y'y'^^biy) = o. (5) 

Furthermore, in RNC the Taylor coefficients of the metric tensor can be found in terms of the covariant 
derivatives of the Riemann curvature tensor |221 EH 125] as 



9ab (y) = Sab + ^Racdb (0) yV + 7:^ eRacdb (0) y'y'y" + -Racdf (0) <,, (0) yVy'y'' 
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^Racdb (0) yY + ^^eRacdb (0) l/"l/V + -^-^acdf y^ ) -t ^^^ 



20 



where y denotes ^{q), the RNC, and q belong to the same patch of p. See appendix A for a derivation of 
this series expansion. Using these coordinates, the geodesic curves look like straight lines passing through 
the point p. Indeed, using the gauge condition ([5]) and the geodesic equation it is easy to figure that out. 
Therefore, the geodesic curve can be written as y"" (s) = ^"s, where s is the geodesic distance and ^" are 



^rV = W''' {db9ec + dcgeb - degbc) 



^ bc^ '^ be ^ c6' anO ^ bed ~ C'ci db "d^ cb^ ^ db'^ ee ^ eb'^ de' 



constants f^. Furthermore, as the geodesic curve is parametrized by the arc-lenght, the coeffitients ^" 
satisfy gab^"'^'' = 1, so the square geodesic distance is given by 

s' = gabV^y' = 6aby^y\ (6) 

where the last equality comes from the conditions of the RNC (|4]). This equation is remarkable, because 
the geodesic distance has the same form as in the euclidean geometry. 

3 Diffusion and geometry 

Here, we introduce the simplest model for the brownian motion of a free particle, which takes place on a 
d-dimensional riemannian geometry. This is a direct generalization of the brownian motion on euclidean 
spaces, which basically consists on replacing the euclidean laplacian by the Laplace-Beltrami operator in 
the diffusion equation [91 [lU] 

dP{x,x',t) nA pr ' +\ 
g-^ = DAgP{x,x ,t), 

P{x,x',0) = —6'^'^^ix-x'). (7) 

Here, P {x,x',t) dv is the probability to find the diffusing particle in the volume element dv = y/gd'^x 
when the particle started to move at x'. The probability distribution P{x,x',t) is normalized with respect 
to the volume v of the manifold and D is the diffusion coefficient. The operator A^ is the Laplace-Beltrami 
operator, which is defined by 

\- = -^d, iVag^'db ■ ), (8) 

where g = det (gab)- We should point out that when the manifold is not compact, we will require 
for the probability and all its partial derivatives to vanish at the boundary. The formal solution of the 
diffusion equation ([T]) on curved manifolds (see [28, 24j) is given in terms of the Minakshisundaram-Pleijel 
coeffients, which depend on both x and x' . This solution has already been used in order to describe the 
concentration of a diffused substance over curved manifolds in the limit when x — )■ x' [12J. 

Once we have a probability distribution, we want to look at the mean-values of physical observables 
(for example, the mean-square displacement) in order to get information of the brownian motion. For 
scalar functions Vl (defined on the manifold), the expectation values are defined in the standard fashion 

{n{x))^= f dvn{x)P{x,x',t). (9) 

Jm 

Note that {Vl (x))^ also depends on the initial point x' . In principle, it is posible to evaluate the expectation 
values using the formal solution mentioned above. However, using this procedure it may be very involve 
because the Minakshisundaram-Pleijel coeffients, as far as we know, are not known for points x different 
from x'. Here, we use a different strategy, which will be applied only for physical observables well behaved 
under actions of A„. 



3.1 General remark on the short time asymptotics 



Clearly, when n(x) is well behaved under any number of actions of the Laplace-Beltrami operator, its 
expectation value can be expanded in Taylor series in the variable t. The k — th derivative of (fi {x))^ at 
t = can be computed as follows. First, let us compute the first derivative using the diffusion equation. 



dt 



D dv Q{x)AgP{x,x',t) 

Jm 

D / dv Agn{x)P{x,x',t)+ / dvVaJ", 

Jm Jm 



(10) 



where J° is a boundary current given by J° = g°''^db{VLP). Using this procedure it is posible to compute 
the k — th Taylor coeffitient by 



d^{n{x)), 



dt^ 



D^A^n (x) 



(11) 



t=o 



here we dropped all the boundary terms because the probability and its derivatives vanish there. The 
expectation value for our physical observable is then given by the formal series [29] 



{nix)), 
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1! 



2! 



Q(x] 



(12) 



where "H = —DAg. This expression is very useful to access the short time regime of the brownian motion. 
Indeed, given a physical observable which is well behaved under the actions of the Laplace-Beltrami 



operator, its mean-value at the short-time regime can be calculated using (12). In particular, for the plane 



case W^ we want to know the mean-square displacement |x| (where x G M^nd we have chosen x' = 0) 



In this case, the Laplace-Beltrami reduces simply to the laplacian <9"9„. Appliying formula (12) we get the 
standard kinematical Einstein relation for the mean-square displacement: ^|x| ) = 2dDt. Observe, we 
can not compute the mean value for the displacement |x| by this method because it is not well behaved 
under the actions of the laplacian at x = 0. 



4 The mean values of s^ and short time asymptotics 

On curved manifold there are several quantities that can be used to described the brownian motion. For 
instance, in [13J the particle position is given in terms of the parametrization of a manifold embedded in the 
ambient space M'^. For this case, the displacement is given by the euclidean norm of the parametrization. 
However, in |T7] the brownian displacement is defined by the geodesic distance. In addition, using the 
Monge parametrization for a surface we can also define a projected displacement (TU]- Here, nevertheless, 
we will not discern between these quantities, instead, we will stress the fact that all of them represent 
different manifestations of the same phenomenon. The analysis between these observables is beyond the 
scope of this work. 

In this paper, we use the geodesic distance as the definition of the displacement of the particle. As in 
the plane case, this quantity is rotational and traslational invariant. Furthermore, the geodesic distance is 



invariant under general coordinate transformations. Therefore the mean-value of s^ will be expected to be 
invariant under a general coordinate transformation. In what follows we compute the expectation value 
of the square displacement using the Riemann normal frame centered at the point p = (f^^{0). Therefore 
the particle starts to move at the initial condition and the expectation value of s^ can be computed using 



(121 



{''). 






fe=l 



k\ 



(13) 



where the terms Gk = •^gS^| {k = 1,2,3,...) are purely geometric factors. This factors can be 
computed explicitly using the technology of the RNC. We will compute the first three factors, Gq, Gi and 

Using the definition of the Laplace-Beltrami operator ([s]) and s^ = 5aby""]/' , it is not difficult to show 
that, for every coordinate y on the manifold. 



A,.^ 



2d + y-d{\ogg), 



(14) 



where g is the determinant of the metric and y ■ d = y'^da- Clearly, when we evaluate at y = we get 
Gi = 2d. The factors G2 and G3 can be found using the result of the appendix B, summarized as follows, 
if / is a diferentiable function on the manifold then the Laplace-Beltrami operator acting on / at y = 
is simply by d"-daf\ q- Therefore G2 = 9^(AgS^)| ^ and G-s = d"^ (A^s^) | _ So, we need at least 



the second order Taylor expansion of A^s and A^s . Using equation (14) and the Taylor expansion of 
log (7 (see appendix A), we obtain the second order of AgS^, 



\y=0- 



A,s^ 



2rf--i?„6(0)yV + 0(2/3). 



(15) 



Therefore, the second geometric factor is G2 = — f-Rg, where the Ricci scalar curvature is evaluated at 



0. For the third gometric factor G^, let Ag act on equation (14) 

1 



^Is^ 



{da log g) g^% {y ■ d log g) + [d^g''') d,{yd log g) + g'^'^dad, {y ■ d log g) . 



(16) 



Using the perturbative expansion of the inverse metric g and the one of the logg it is not difficult to 
obtain the second order of A^s^ (see appendix A). Therefore, by a straightforward calculation, we get 



G, 



15 



K Kab 



16 

45 



jyabcd / p I p 

-K [rCdbca + -K, 



dcba ) 



16 



V-V" + -(/""A, 



R. 



ah- 



(17) 



In general, for the Gk factor we need to compute the second order pertubative expression for A^ 
For the expectation value of s^ at the short-time regime we have considered only the values k = 
Hence, this value can be written as 



1,2,3. 



(s^)^ = 2dDt - -Rg {Dty 



1 
3! 



15 



ryab p 
K Kab 



16 
45 
16 



j^abcd (^ 



dbca 



Rdcba) 



V^V" + -g'^'/^g 



R. 



ab 



{Dtf + 



As we anticipate in the introduction, the mean-square geodesic distance (i.e. (s^)J is deviated from 
the planar expression by terms which are 0(d) invariant as well as invariant under a general coordinate 
transformation. Clearly, for very short-times, the particle movement is not affected by the curvature of 
the manifold and the standard mean-square displacement is recovered [IS]. This follows from the very 
nature of the manifold (every local neighbourhood looks like a piece of an euclidean space). However, as 
the particle explores away from the local boundaries, the movement is affected by the curvature of the 
manifold. Additionally, the curvature corrections to the planar expression are isometric. Therefore, every 
manifold which is isometric to an euclidean space, will have null curvature effects for this observable. In 
particular, this is the case for the developable surfaces [17j . 

4.1 Example: the brownian motion on the n-sphere 



Here, the idea is to perform a cross-check calculation in order to validate equation (18). In particular, we 
compute the mean-square geodesic distance at the short-time regime, when the brownian motion takes 
place on a n-dimensional sphere. The hypersphere S" of radius R is defined by 



5" 



{xG 



r+i : x^ 



R^}. 



(19) 



This manifold can be parametrized using the local coordinates (9o,9i,--- ,9„-i), where 9q takes values in 
[0, tt] whereas the remainder coordinates in [0,27r). In this example, we use a parametrization concerning 
to the ambient space ]R"+^ by using the functions x = (a;o, ■ ■ ■ ,x„) given by 



xq = R cos 6*0 

Xi = Rsm9QCOs9i 



J^r, 



R sin 9q sin 9i- ■ ■ sin 9n-i- 



(20) 



Clearly, x = i?n, where n is the unit normal vector pointing outward from hypersurface of S^. The metric 
tensor can be computed using gab = da^ ■ dbX and it can be written in a matrix form as 



/I 



9ab 







\ 



sin^ 9n 



sin^ 9o sin^ 9i 



\ 







(21) 



sin^ 6^0 ■■ ■ sin^ 9n-i j 



The square root of the metric tensor determinant is y/g = R" sin"^^ 9q sin"~^ 9i- ■ ■ sin6'„_i. It is remar- 
kable that using extrinsic geometry we can easily compute the Riemann curvature tensor. Indeed, using 
the Gauss-Codazzi equations, Rabcd = KacKbd — KadKbc, where Kab = t^aX ■ db^v are the components of 
the second fundamental form [20], the Riemann curvature tensor is 



R. 



abed 



R" 



{gacdbd — Qaddbc 



(22) 



The Ricci tensor and the scalar curvature are then given by 



R, 



ab 



n — I 



-Qah, 



Rn 



n{n — 1) 
R^ ' 



(23) 



respectively. Substituting equations (22) and (23) into the expressions of Gi, G2, and G3, obtained above, 
we find 



n -o^ n - ^ ^(^~^) n -'^ n(n-l)(n-3) 



(24) 



Now, in order to perform an indepedient calculation we compute the mean-square geodesic distance using 



the equation (12). For practical purpose, we use the geodesic curve starting at the north pole defined as the 



point X = (1, 0, ■ ■ ■ , 0) of 5". It is not difficult to show that Qq = j^s and 9j = (with j = 1,- ■ ■ ,n — l) 
is this geodesic curve. The geodesic distance is simply given by s = ROq. 
In this case, the Laplace-Beltrami operator on S" can be splitted as 



n — 1 



A 



^ R^ sin"-i ft 






(25) 



where 5o = ^ and 9 



Clearly, the actions of the Laplace-Beltrami operator on s^ will involve only 



the first term of (25). The first, second and third action of A^ on s^ are given by 

A^s^ = 2(n- 1)^0 cot ^0 + 2, 

^2 [n - 1 + (n - 3)(1 + cot^ ^0) (^0 cot Oq - 1) 



AV = -2(--l) 



3„2 



^> 



2(n-l)(n-3) 



{{n-1) cot ^0(1 + cot^ 60) [cot ^0(3 - 00 cot Oq) - 0o{l + cot^ 6^ 
+ 4(1 + cot^ ^0) [(^0 cot 00 - 1)(2 cot^ ^0 + 1) - cot ^o(cot Oq - 9o{l + cot^ ^0))] } • 



i?4 



(26) 



Therefore taking the limit when ^0 — ^ we get (24). The mean-square displacement is then given by 



[s')^ = 2nDt - 



2n{n-l) 2 4: n{n-l){n-3) 



i?2 



(Dty + 



45 



i?4 



{Dtf + ■ ■ 



(27) 



This equation is the desired result (18) for the particular case of the n-sphere. 



4.2 The constant curvature spaces 

For the constant curvature manifolds, the Riemann curvature tensor is given by [30] 



R 



■abed 



d{d-l] 



cgbd — gadgbc), 



(28) 



where Rg is a constant that can be either positive, negative or zero. It is known [3U] that the only three 
solutions for constant curvature are the d-dimensional sphere when Rg > 0, the d-dimensional hyperboloid 
when Rg < and the euclidean space for Rg = 0. For these cases, the mean-square geodesic distance is 
given by 



s^)^ = 2dDt - ^Rg {Dtf + ^ 



rf-3 



45 d{d - 1 



< {Dtf + ... 



(29) 



For theses cases, the curvature effects depends on the sign of the curvature Rg. For the hyperbolic spaces, 
the MSD deviates from the planar result by an increasing monotonic function (in time) whereas for the 
spherical space by a decreasing monotonic function. In other words, the geometry of the space affects the 
brownian motion in such a way that the particle's diffusion is accelerated for i?^ < or decelerated for 
Rg > 0; see figure (jTl). These same results were also observed in the two-dimensional cases in [18J. 




Figure 1: The time evolution of the MSD for the spherical Rg > 0, hyperbolic Rg < and euclidean spaces 



Rg=0. 



4.3 The two-dimensional case 

The two-dimensional case is the most relevant one for the diffusion on biological membranes. For this 
case, the Riemann curvature tensor is given by 



R, 



abed 



K {gacQbd - 9ad9bc 



(30) 



where K = Rg/2 is the Gaussian curvature of the surface [21]. Note that K is not necessarily a constant. 
The mean-square geodesic distance, at the short-time regime, for an arbitrary surface of Gaussian curvature 
K is given by 



[s')^ = 4Dt--K{Dty 



15 



-K^ + 2A„K 



{Dtf + 



(31) 




Figure 2: A Torus divided by three regions, where K > , K < (i.e the area gridded) and K 
(parallels circular curves where = 7r/2 and 6 = 37r/2). 



In what follows, we derive few consequences of last equation when the diffusion takes place on a Torus. 
For this case, the metric is given by ds'^ = r'^dO'^ + {a + r cos 6)^ d^p^ (with < 6^, </) < 27r), where a and 
r are the two radii. The Gaussian curvature is given by 

K= ^ '''' ^^ . (32) 

r[a + r cost/) 

At figure ([2]) three regions on the Torus are shown. These regions are defined by the conditions K > 0, 
K < and K = 0, respectivel^R It is shown at figure dsl) how the time evolution of the mean-square 
displacement depends on the initial position of the particle on the Torus. Clearly, the particle's diffusion 
decelerates faster at the region where K < than the way how the particle's diffusion proceeds at the 
region where K > 0. On the contrary, at the parallels [K = 0) the particle's diffusion accelerates; in fact, 
the MSD in the parallels is given by 

Ifi 
(A = 4^^+I^^(^^)'- (33) 

In general, the biological membranes have a wide range of morphologies; from discoidal shapes to 
catenoidal ones [31]. Also, we can consider wavy surface like membranes with microvilli [9J, elliptic 
paraboloid, hyperbolic paraboloid [18J and periodic nodal surfaces [13J, where the diffusion is clearly affec- 
ted by the geometry. The curvature effects on the diffusion on these surface can be quantified using 



equation (31) 



"^ Clearly, the Gaussian curvature vanishes along the parallels {6 = tt/2 and 9 ~ 3tt/2). In the region given by 
7r/2 < 6 < 37r/2, K is negative whereas in the one given hy < 9 < 7r/2 or tt/2 < 9 < 37r/2, K is positive. 
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Figure 3: The time evolution of (s^)^ versus t and 9, where the area gridded represents the region where 
K <0. 



5 Conclusions and perspectives 

In this paper, we have studied the brownian motion over a d-dimensional curved manifold. The geodesic 
distance is used as the displacement of a particle diffusing on the space. Using a Riemann normal frame we 
have derived a general formula for the mean-square displacement at the short-time regime. This formula 
reproduces the standard result for very short-times [16J and it is given by 0(d) isometric covariant terms 
depending upon the Riemann curvature tensor. In particular, we have explored the diffusion over constant 
curvature manifolds, where we have shown that diffusion accelerates for the hyperbolic case, whereas it 
decelerates for the spherical one. In addition, we have discussed the two dimensional case, as it is relevant 
for the diffusion on biological membranes. In particular, the behaviour of the diffusion on a Torus can be 
classified according to the region (defined hy K < 0, K > or K = 0) where the particle starts to move. 

As we have mentioned, there are several physical observables that describe the motion of a particle 
diffusing on a curved space. These observables quantify the curvature effects in different ways. For 
instance, as it was discussed at [IT] the mean-square geodesic distance has null curvature effects for 
the diffusion on developable surfaces, whereas using the parametrization displacement [13J it is clear that 
there are curvature effects. Indeed, different physical observables have different manifestations of the same 
phenomenon. It is then a rather natural question to ask what is the geometrical and physical content of 
these physical observables for the brownian motion on curved manifolds. 

In a different direction, we can explore the large-time regime for the mean-values for the case of 
compact manifolds. In this particular case, the expectation values will have a bounded above limit as a 
consequence of the compactness of the manifold [29j. 
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7 Appendix A. Riemann normal coordinates 

7.1 The gauge condition 

As it is mentioned at [25] the second condition of Q is equivalent to the gauge condition (^. In order 
to prove this we use the Christoffel symbols, therefore we have 

y'^y'T^.iy) = ^^- {2y'^y''daget - y^y^gab) ■ (34) 

Now, the first derivative on the second condition of Q is gab + y^dageb = ^ab, then y'^y^dageb = 0. 



Similarly, the second term of (34) vanishes 



7.2 An afRne connection expansion 

In order to derive (|6]) we follow the same procedure presented at |22]- Let w be a matrix one-form 
connection given by a; = ojadx"', where the matrix one-form component {uJa)\ = ^^ac '^ given by the 
Christoffel symbols. Let us consider the Lie derivative Cx along the radial vector X = y ■ d acting on cu, 

Cxto =[{l + y d)LOa] dy'', (35) 

where y"- are the Riemann normal coordinates. The Lie derivative can be written as 

Cx = dix + ixd, (36) 

where d and ix are the exterior and interior derivative, respectively [27j. Observe, that the gauge condition 
(^ can be written as ix(^ = 0. Therefore using the Cartan structure equation for the curvature, R = 

dcu + cu ALU, is easy to get 

Cx<-^ = 'ixR — ixiy /\ ^) 
= ixR 
= y^TZabdy". (37) 



Equating (35) and (37), we get the condition {1 + y ■ d)ua = y"''R-ab, where {TZabYd = ^'^dab '^ the 



Riemann curvature tensor . A Taylor expansion centered at in both sides of this condition gives 

{y ■ df 

k\{k + 2) 



^''bc = Yl TmrT^y''^''cdb, (38) 



k=0 

where Riemann curvature tensor is evaluated at 
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7.3 A vielbein expansion 



As it is point out in |22], it is sufficient to compute the Taylor expansion of the vielbein 9 = Oadx"" in 
order to find ([6]). 6 satisfy the free torsion condition dO + uj /\6 = {) and the metric can be re-written in 
terms of the vielbein as gab = 0\0\Sij. Using (36) and the gauge condition, the Lie derivative on 9 is 
given by 

Cx9 = {ix9)uJ + dixO. (39) 

As 6 is one-form, ixO is scalar. Therefore the second Lie derivative of last expression is 

CxCxe = X [1x9] uj + [1x9) Cxuj + CxdixO. (40) 

The additional gauge condition y°-gab {y) = V^^ab is equivalent to choose y°'9l (y) = y"'5l, then X [ix9] = 



Now using equation (37), we find 



ix6 and Cxdixd = dixO. Combining equation (39) with (40), we get Cx{Cx — I) 9 = {ixO) Cx^ 



Cx{C 



X 



(ix9)ixR- 



(41) 



The condition for the vielbein is then given by {y ■ d + 1) (y ■ 3)9^ = y'-' y^ R!' jijj)^ ■ By Taylor expansion 
centered at in both sides of this condition we get 



ei 



yde\ 



s. 



^ {y ■ df-' [R\Jl\ 



--/l.fc(fc + l)(fc_2)! 







(42) 



where i?*..^ = y°-y R\^. This system of equations can be solved iteratively and they give the following 
expression for the vielbein until fourth order in y°' 

o\ = ^1 + \R\.a + l^y ■ dR\.a + ^{y df R\.a + ^R'-cR^'-a + ■■■ (43) 

The metric (pj) is then obtained using gab = d\Q\^ij- 



7.4 Determinant of the metric 

A Taylor expansion of the determinant of the metric g = det gab can be obtained using the identity 
logdet gab = tr\oggab- Since the metric ^ can be written as gab = ^ab + ^ab (y), its logarithm is given 
by log(l + A) ^ A — |A^. The trace of this term is then given by 

-\ogg = l^Rab (0) y'^y' + l^^Rab (0) y^yV + ^R'cdj (0) i?,^,, (0) yVyV 



6 

+ ^VeV/i?,,(0)l/Vz/V + 



(44) 
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7.5 Geometric G3 factor 

In order to get the geometric factor 6*3 = A^s^, we need at least the Taylor expansion of A^s^ until 
second order 0{y'^). By straighforward calculation this expansion is given by 



Ajs^ = -\Rg- y^'daRg - 2y^V'Rab + \R:RadyV + \ R^,' Rai>yV - ^ {-2R}Rf^,yy'' 

+ RZfKhavV + R'^fR^ .Hav'y' + R\fR\ayV + Rac '^Rf,H ""fy'} - 1 {v'Rca) y'y 



4„ „_ __.„ 4„„„ _ 4„„ _ _ 4^_^^ 

2 

- ^(VeV/i?,)i/V-^(V^V/i?ed)l/V + 0(i/='). (45) 



8 Appendix B. A^ on a scalar function 

Let f] : M — 7> M be a scalar differentiable function on the manifold. Let us take a Riemann normal system 
of coordinates centered at y"- = 0, then the Laplace-Beltrami operator acting on Vl is given by the second 
derivative of ^2, i.e., 

A„n| n=^^^l n- (46) 

y ij/=o ij/=o ^ ' 

In order to show this result we split the inverse metric and the square root of the metric determinant as 

gab ^ 5^^ + K^\y) 

V9 = 1 + A(y), (47) 

where A"* (y) and A (y) can be found in terms of covariant derivative of the Riemann curvature as we 
have seen above. This functions satisfy the following properties: A"^ (0) = 0, (9c A"* (0) = 0, A (0) = 0, 
and dc\ (0) = 0. Thus the action of A^ on fi can be written as 

AgQ = {d^Q + daX{y)d''Q + X{y)d^Q + daA''''{y)dbQ + A''''{y)dadbn 

+ daX (y) A"' (y) d,n + A {y) daA'^' (y) d,n + A (y) A'^' (y) dad,n}l^^ 

(48) 

Then after the substitution y = 0, it is obtained the desired result. 
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